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Abstract 



Fixed point actions for free and interacting staggered lattice fermions 
are constructed by iterating renormalization group transformations. At 
large N the fixed point action for the Gross-Neveu model is a perfect 
action in the sense of Hasenfratz and Niedermayer, i.e. cut-off effects 
are completely eliminated. In particular, the fermionic 1-particle energy 
spectrum of the lattice theory is identical with the one of the continuum 
even for arbitrarily small correlation lengths. The cut-off effects of the 
chiral condensate are eliminated using a perfect operator. 
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1 Introduction 



Cut-off effects are the main source of systematic errors in numerical simulations 
of lattice field theories. Artifacts due to a finite lattice spacing a vanish in the 
continuum limit when the correlation length diverges in lattice units. In bosonic 
theories the lattice artifacts are usually of 0(a 2 ) and in fermionic theories they 
are of 0(a). Hence they go to zero rather slowly as the continuum limit is ap- 
proached. In practice it is very difficult to work at large correlation lengths mostly 
because of critical slowing down. When Wilson introduced lattice field theory his 
idea was to use the renormalization group to map the critical continuum theory 
— defined at a fixed point of the renormalization group — to a noncritical theory 
with small correlation length which could then be solved numerically [I]. This re- 
quires to approach the renormalized trajectory of theories connected to the fixed 
point by an infinite number of exact renormalization group steps. Theories on this 
trajectory have perfect lattice actions — their spectrum is completely free of cut-off 
effects. Unfortunately, attempts to locate the renormalized trajectory using real 
space renormalization group techniques have not been successful. Therefore numer- 
ical simulations must approach the continuum limit the hard way, going to larger 
and larger correlation lengths until scaling is finally observed. 

The idea of Symanzik's improvement program [2] is to systematically approach 
the renormalized trajectory by adding irrelevant operators of a given dimension to 
the standard lattice action. Their coefficients are fixed by demanding that physical 
quantities become cut-off independent up to a given order in a. As the desired order 
of improvement is increased, more and more operators of higher and higher dimen- 
sion must be taken into account. In asymptotically free theories the coefficients of the 
various operators can be computed in weak coupling perturbation theory. The lead- 
ing contribution defines classical (tree level) improvement, while the higher orders 
represent quantum (re-loop) improvement. For gauge theories Liischer and Weisz 
[3] have introduced the concept of on-shell improvement and they have constructed 
the tree level improved action for Yang-Mills theory to 0(a 2 ). Sheikholeslami and 
Wohlert [4] have applied similar ideas to QCD and their 0(a) tree level improved (so 
called clover) action is now often used in numerical simulations. In the Gross-Neveu 
model Wetzel [5] has implemented Symanzik's improvement program in the large N 
limit. 

Recently, Hasenfratz and Niedermayer [6] have taken a more radical approach 
to eliminate cut-off effects. They realized that perfect actions can be constructed 
explicitly for asymptotically free field theories. Their method is nonperturbative 
and it results in the full fixed point action by working with an, in principle, infinite 
number of irrelevant operators. In practice one is of course limited to a truncated 
(but still large) set of operators. The fixed point action is a perfect classical action 
and it should resemble a tree level Symanzik improved action to all orders in a. 
Hasenfratz and Niedermayer also discuss how to derive perfect quantum actions 
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which are located on the renormalized trajectory (and which should correspond to 
all-loop improvement to all orders in a). However, in practice this seems not to be 
really necessary, at least in the 2-d 0(3) nonlinear cr-model that they investigated 
numerically. Using the perfect classical action at correlation lengths as small as 5a 
the cut-off effects were smaller than the statistical errors and were hence practically 
eliminated. At present one does not really understand why the perfect classical 
action works so well at the quantum level even at these extremely small correlation 
lengths. 

The method of ref. [6] applies to any asymptotically free theory. Of course, we are 
most interested in QCD, and it is a great challenge to construct its perfect action. 
However, this needs some preparation and cannot be done in one step. A natural 
first step is the construction of a perfect action for pure Yang-Mills theory which is in 
progress [7]. In the second step one would include the quarks. We think that it may 
be helpful to gain experience with fixed point actions for fermions also in simpler 
settings. The simplest asymptotically free theory with fermions is the Gross-Neveu 
model in two dimensions [8]. In the continuum the model with N flavors is defined 
by the euclidean action 

. N ri N 

8 = 1 8 = 1 

For later convenience we rescale the fields to 

? = VGV, X l = ^\ (1-2) 

and we introduce a real valued auxiliary field $ to linearize the 4- Fermi interaction 
by a Yukawa coupling. Then the action turns into 

1 l f N . .1 N . . 

7^[x, x^ = q] d " x ill r-rM + + E (1-3) 

8 = 1 8 = 1 

and the 4- Fermi coupling constant appears as a global prefactor. The model has an 
0(2N) (x) Z(2) chiral symmetry that gets spontaneously broken to 0(2N) resulting 
in a dynamically generated fermion mass mj [8]. The model is asymptotically free, 
i.e. in a lattice formulation the continuum limit corresponds to G — > 0, and for 
small values of G one expects 

amj(G) ~ exp(—TT /G(N — I)). (1.4) 

Note that the model with a single flavor (N = 1) is identical with the Thirring 
model which has a vanishing ^-function. Also the N = 2 model is special. It has 
been conjectured that it is equivalent to two decoupled sine- Gordon models. The 
large N limit is taken such that g = GN is kept fixed. Then asymptotic scaling 
corresponds to 

am f (g) ~ exp(-wfg). (1.5) 
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The lattice Gross-Neveu model has been studied e.g. in refs.[9, 10, 11] using stag- 
gered fermions. Then the number of flavors is even. On the lattice the 0(2N) ® Z(2) 
symmetry of the continuum theory is explicitly broken to U(N/2) e=0 Cg) Z(2). It is 
essential that the Z(2) symmetry remains intact, because then the dynamics of its 
spontaneous breakdown is not obscured by lattice effects. 

Here we investigate fixed point actions for the Gross-Neveu model in the large 
N limit, where the whole calculation can be done analytically. It turns out that at 
N = oo the fixed point action — which is the perfect classical action — is in fact also 
a perfect quantum action. In particular, the energy spectrum of the lattice theory is 
identical with the one of the continuum even for arbitrarily small correlation lengths. 
This shows that the success of the method of Hasenfratz and Niedermayer is not 
limited to the (9(3)-model, and it hopefully indicates that it will also work for QCD. 
We also investigate a small field approximation to the fixed point action which turns 
out not to be a perfect action. This underlines that the nonperturbative approach 
of ref. [6] is really necessary, and that one should not seek short cuts using small field 
approximations, although this is usually much simpler than determining the full fixed 
point action. Moreover, we construct a perfect operator for the chiral condensate 
to demonstrate that its cut-off effects can also be eliminated. However, this does 
not happen automatically just by using a perfect action. Then we demonstrate that 
with the perfect action asymptotic scaling sets in earlier than with the standard 
action. Still, for the perfect action asymptotic scaling is not perfect. Asymptotic 
scaling is not really what one should ask for. It is an unphysical issue because it 
involves the bare coupling defined at the cut-off scale. The real issue is scaling of 
ratios of physical quantities — and scaling is perfect for a perfect action. Finally, 
we investigate the renormalized trajectory in order to understand why the perfect 
classical (fixed point) action is even a perfect quantum action. It turns out that 
this is the case only at N = oo. For finite N one expects cut-off effects at finite 
correlation lengths even with the perfect classical action. Still, the fact that it 
becomes perfect at large N may explain why the fixed point action was practically 
perfect in the 0(3) model study of Hasenfratz and Niedermayer. Of course, this 
argument assumes that the model at N = 3 does not deviate significantly from the 
large N limit. 

In section 2 we investigate the fixed point actions for free staggered fermions 
using a renormalization group transformation suggested by Kalkreuter, Mack and 
Speh [12] that is consistent with the staggered symmetries. We generalize their 
renormalization group transformation in a way that allows us to optimize the lo- 
cality of the fixed point action. In section 3 we turn to the 2-d Gross-Neveu model 
by switching on a 4-Fermi coupling which is linearized by a Yukawa coupling to an 
auxiliary scalar field. Blocking both the fermion and the auxiliary scalar field we 
determine the fixed point action of the coupled system in a small field approxima- 
tion. Section 4 deals with the large N limit. There we can go beyond the small 
field approximation for the zero mode of the auxiliary scalar field, which allows us 
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to determine the perfect classical action analytically. We also construct a perfect 
operator for the chiral condensate. In section 5 the Gross-Neveu model is solved 
in the large N limit and it is shown that the perfect classical action is in fact also 
a perfect quantum action. In particular, the cut-off effects of the energy spectrum 
and of the chiral condensate are completely eliminated. We also investigate the 
renormalized trajectory to answer the question why the fixed point action is even a 
perfect quantum action. Finally, we draw some conclusions in section 6. 



2 Fixed point actions for free staggered fermions 

The results of this section have been presented before by two of the authors [13]. 
A similar study for Wilson fermions was performed in ref. [14]. In the Gross-Neveu 
model chiral symmetry is spontaneously broken. Since Wilson fermions break chiral 
symmetry explicitly it is easier to investigate questions of spontaneous chiral sym- 
metry breaking using staggered fermions. They are protected against perturbative 
radiative mass corrections by a remnant of chiral symmetry. Staggered fermions 
are also easier to use in numerical simulations because the location of their critical 
surface is determined by chiral symmetry. 

In two dimensions staggered fermions represent two flavors of Dirac fermions 
in the continuum limit. Staggered fermions can be described by one-component 
Grassmann variables living on a quadratic lattice with spacing a/2. Here we label 
the staggered fermion variables x l x i X x by t wo indices x and z, where x are the centers 
of 2 X 2 disjoint blocks and i £ {1,2,3,4} is a pseudoflavor index that determines 
the position within a block. The pseudoflavors 1 and 2 correspond to the lower left 
and right corners while the pseudoflavors 3 and 4 represent the upper left and right 
corners of the block. The block centers x form a quadratic lattice with spacing a. 
For technical reasons we locate the fermions at the 2x2 block centers x and not at 
the sites of the lattice with spacing a/2. From now on we set a = 1. In general the 
action of the free theory may be parametrized as 



Here p(z) is a 4 X 4 matrix in pseudoflavor space that depends on the difference 
vector z = x — y between two block centers. Staggered fermions have various 
symmetries, among them a U(l) e <E> U(l) remnant of chiral invariance, an analog 
of charge conjugation, and pseudoflavor transformations that correspond to shift 
symmetries on the lattice with spacing 1/2. The symmetries impose constraints on 
the matrix p(z). The U(l) e (£)U(l) symmetry implies pu(z) = 0, pi 4 (z) = P4i(z) = 0, 
P2?>{z) = P32(z) = 0, charge conjugation leads to pi^—z) = —p 1 i(z) } and the shift 
symmetries give in addition P24(z) = — p\3(z) } P3 4 (z) = pi 2 (z). This leads to the 



S[X, X] = Y,Y1 XxPij(x - y)Xy = Xxp{x - y)xy 
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structure 



( Pl (z) p 2 (z) \ 

-Pi(-z) -p 2 (z) 

-p 2 (-z) Pl (z) 

V p 2 (-z) - Pl (-z) J 



(2.2) 



The shift symmetries restrict the matrix elements even further by 

P„{-z) = -p„{z + p), (2.3) 

where p is the unit vector in //-direction. The standard action is characterized by 
p(z) = p s (z) with 

^(0) = -l, Pl(p) = l. (2.4) 

All the other values are zero. The actions that we will generate by renormalization 
group transformations are more general but they are consistent with the symmetry 
requirements summarized in eqs. (2.2,2.3). Note that we have not introduced a chiral 
symmetry breaking fermion mass term. The model represents massless free fermions 
with infinite correlation length and hence it corresponds to a critical lattice field 
theory. 

The standard action represents a particularly simple point on the critical surface 
but, except for this, not a very interesting one. Here we are interested in fixed 
points of the renormalization group because this is were the continuum limit of a 
lattice field theory is defined. To locate the fixed points we start at any point on 
the critical surface — in this case at the point defined by the standard action — 
and we iterate renormalization group transformations. The location of a fixed point 
will in general depend on the renormalization group transformation that we choose, 
and we will use this freedom to optimize the fixed point action's locality. To be 
consistent with the symmetries of the problem it is important that the renormaliza- 
tion group transformation respects the pseudoflavor structure of staggered fermions. 
Kalkreuter, Mack and Speh have proposed a suitable blocking scheme with blocking 
factor 3 [12]. To respect the staggered fermion symmetries the blocking factor must 
be odd. The blocked coarse lattice has spacing 3. Hence nine 2x2 block centers x 
form a block that is associated with a new 2x2 block center x' on the coarse lattice. 
We denote this by x £ x 1 . Each pseudoflavor is blocked individually such that each 
\ l x on the original fine lattice contributes to exactly one x'x' on the blocked coarse 
lattice. The block transformation is illustrated in fig.l. First we apply a ^-function 
renormalization group transformation such that after one blocking step the effective 
action S'[x\x'} ls given by 

l „ r , l 



1 t t r 1 

exp(-— S"[x',x']) = J VxVxexp(-— S[x,x\) 



n$(x'x> xDHxx* ~^J2xl 

x ,i x£x x£x 
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Figure 1: Blocking of a 2-d lattice consistent with the staggered symmetries. The 
symbols •, o, o, * represent the pseudoflavors 1, 2, 3, 4 on the lattice with spacing 
1/2. Four 3x3 blocks of the various pseudoflavors are also shown. Their centers 
forming the blocked lattice are marked with a large circle. 

= j VxDxVrjVri exp(-^S[x, x]) 

x ex P {££[(x£ - A j2 xi)vi> + vl>(x':> - 4 E xi)]}- 

i O _ / O - i 

x % x£x x£x 

(2.5) 

We have already introduced the 4-Fermi coupling G of the Gross-Neveu model. 
There due to asymptotic freedom the continuum limit (and hence the critical surface) 
is at G — > 0. Then, as Hasenfratz and Niedermayer pointed out, finding the fixed 
point action is a classical saddle point problem. For each pseudoflavor we have 
introduced a coarse lattice auxiliary Grassmann field ff xll rf x ,. The parameter b 
renormalizes the fermion field. It must be fixed appropriately in order to reach a fixed 
point of the renormalization group. It is important to note that the renormalization 
group transformation leaves the partition function — and hence the physics — 
invariant. This follows immediately when one integrates out the blocked variables 

Z' = j Vx'V X 'exp(-^S'[x',x']) = j V X V X eM-^S%x]) = Z. (2.6) 

Before we proceed with the blocking step we generalize the renormalization group 
transformation by introducing a kinetic term for the auxiliary Grassmann field 

1 f 1 

exp(-— S"[x',x']) = J VxT>xT>T]VT]exp(-— S[x,x\) 

x' % x£x' x£x' 
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x eM--Ev x >p s (^—T^)Vy>)- (2-7) 
c 4 



Here p s (z) is the matrix defined in eq.(2.4) that characterizes the standard action. 
The parameter c will be used later to optimize the locality of the fixed point action. 
It is easy to show that with the extra kinetic term for the fermionic auxiliary field the 
renormalization group transformation still leaves the partition function invariant. It 
should be noted that this is not the case for gauge theories. Hence, there one is 
limited to c = oo. Of course, one could also introduce a mass term for the auxiliary 
fermionic field. This would, however, break the remnant chiral symmetry and would 
therefore contradict our strategy. This is in contrast to Wilson fermions where chiral 
symmetry is completely broken already by the action. In that case it is natural to 
introduce additional chiral breaking via the renormalization group transformation 
[14]. 

To perform the renormalization group step it is useful to go to momentum space. 
Then all variables are replaced by their Fourier transforms, e.g. 

P(z) = f B d 2 k p(k) exp(ikz). (2.8) 

Here we integrate over the Brillouin zone B =] — 7r,7r] 2 . The quantities defined on 
the blocked lattice are integrated over a smaller Brillouin zone B' =] — 7r/3,7r/3] 2 . 
In momentum space eq.(2.7) takes the form 



exp(- — S"[x',x']) = j VxVxVrjVriexpi- — J^y J B d 2 k x(-k)p(k)x(k) 
+ {^j j B> d 2 k W{-k)r,{k) + fj(-W(k)] 

~ i d 2 kr](-ky(M) V (k)}. (2.9) 

c \ZttJ Jb' 



We have defined 

m = m-^j:ii(k + 2 -f)D(k + 2 -f)x(k+ 2 -f) } 

n-k) = ^- k )-±^ x (-k-^)D(-k-^)Il(k + ^-). (2.10) 

The sums extend over integer vectors with components £ {1,2,3}. We have also 
introduced a matrix D(k) diagonal in pseudoflavor space that is given by 

% % 
Du(k) = exp & 2 ), D 22 {k) = exp -(k x - k 2 ), 

D 33 (k) = exp'-i-h + h), D 44 (k) = exp t -(k 1 + k 2 ). (2.11) 
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It arises because in our notation the fermionic variables are located at the 2x2 
block centers x forming the lattice with spacing 1, and not at the sites of the lattice 
with spacing 1/2. Finally, the 3x3 block average gives rise to the function H(k) = 
EL |(l + 2 cos k v ). To obtain the effective action we first integrate out the fermionic 
auxiliary field 

1 f 1 1 f 

exp(-— S"[x ,x']) = J ^X^Xew{-^J^y j B d 2 kx(-k)p(k)x(k) 

+ ^(^) 2 J B /k^\-k)a'(3k)^(k)}, 



:2.i2) 



where a s (k) = p s (k)~ 1 . Next we solve the saddle point problem, i.e. we derive 
classical equations of motion for the fermionic variables by varying the exponent 
in eq.(2.12) with respect to x{ — k). Keeping in mind that if)'(-k) also depends on 
x{ — k) one finds the classical solution 

Xc (k) = -bca(k)D(-k)Tl(k)a s (3k)ip'(k) } (2.13) 

where a(k) = p(k)~ 1 . Multiplying eq.(2.13) by D(k)Tl(k) and summing it over the 
momenta k + 2ttI/3 with £ {1, 2, 3} one obtains 

tjj'(k) = x'(k) + cu}(k)a s (3k)^'(k), (2.14) 

with the matrix 

/,x b 2 27r/. ., 27r/. „. , 27r/.„., 27r/, 

<"(*) = ^J2 D ( k + —H k + —W~ k ~ —W k + — ) ■ ( 2 - 15 ) 

^From eq.(2.14) we infer 

4>'(k) = (1 - cuji^a'iUyyWk). (2.16) 

Similar relations hold for x c ( — k) and if)'(-k). Inserting these as well as eq.(2.13) and 
eq.(2.16) back into the exponent of eq.(2.12) one obtains after some manipulations 

S'\x', x] = ~c {£) j Bi d 2 k X '(-k)a s (U)(l - co;(*X(3*))~ Y(*)- (2.17) 
Hence one identifies 

a'^k)- 1 = p'(3k) = -ca s (3k)(l - cto{k)a s {Sk))' 1 . (2.18) 
Using eqs.(2.16,2.18) we rewrite eq.(2.13) as 

Xc (k) = bai^Di-^Ui^a'iUyWk). (2.19) 

Similarly, one derives 

Xc(-k) = bx'i-kyiM^ttikjDikjaik). (2.20) 
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These equations will be useful later because they allow to determine the fixed point 
action for the interacting theory in a simple manner. Now we return to eq.(2.18) 
which we rewrite as a'(k) = — ^p s (k) such that 

<* (*) = 32 IX^g— W — ^— )n(^— ) 2 - - P '(k). (2.2i) 

To eliminate the D-factors we introduce the matrix a(k) = D( — k) 1 / 2 a(k)D(k) 1 / 2 
which has the structure 1 



a(k) 



( a x {k) a 2 (k) \ 

a^k) -a 2 {k) 

a 2 (k) ai(&) 

V -a 2 (k) &!(k) j 



(2.22) 



In coordinate space the multiplication with D(k) 1 / 2 corresponds to a pseudoflavor 
dependent shift of the fermionic variables by ±1/4 of the lattice spacing in each 
direction. This shifts the fermions away from the centers x of the 2x2 blocks (that 
form the lattice with spacing 1) to the sites of the lattice with spacing 1/2 that is 
illustrated in fig. 1 . Generally, all quantities with tildes (there will be more defined 
later) refer to these shifted fermions. Eq.(2.21) turns into the recursion relation 

//7 , b 2 ^ ,& + 27rf w w tt / 2sin(kj2) \ 2 1 f , 
a'(k) = —Y&J— -l' M TT — - — - — — — --— - — - ) +-ik u . 2.23 

We define p s (k) = D(-k) 1 / 2 p s (k)D(k) 1 / 2 and we have used = —ik^ with 

k^ = 2sin(A; |U /2) as well as 1 + 2cos(k l/ ) = sin(3A; J/ /2)/ sin(k l/ /2). The sign factor 
( — 1)' M in eq.(2.23) arises after the D-factors have canceled. Now we iterate eq.(2.23) 
and after n renormalization group steps we find 

» [) \3 2 ) Y 3 " ^ \3-2sin((^ + 2wl v )/&2) J 

l-(b 2 /3 3 ) n f 

with 1^ £ {1,2,..., 3 n }. In the limit n — > oo only the small k behavior of a^k) close 
to its pole in B is important. Here we have a^k) = ik^/k 2 ~ ik^/k 2 such that after 
an infinite number of renormalization group steps 



1 In ref.[13] we forgot to put the tilde on a(k). 
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The limit is finite — i.e. a nontrivial fixed point is reached — only il (b 2 /3 2 ) n 3 n = 1 
which implies b = y/3. This is what one would expect on dimensional grounds. The 
canonical dimension of a fermion field in two dimensions is d x = 1/2. In order to 
renormalize the fermion field by the appropriate amount, b should be given by 3 d *. 
At the fixed point one obtains 

x-^. k,, + 2tvL, , w T-r / 2 sin(&,,/2)\ 2 9 ? , , 

- . e ^frw'- 1 ' " n (x^J + £ (^) 

Now the parameter c is used to optimize the locality of the fixed point action. For 
this purpose we consider field configurations that are constant in the 2-direction, 
i.e. in momentum space they have k 2 = 0. For these effectively 1-dimensional 
configurations it is sufficient to consider «*(&!, 0). One finds o^&ijO) = and 

n • 1 , w /2sin(&i/2)\ 2 9 . ? .A 1 ? 9- 

(2.27) 

For c = 9 we have /?^(&i,0) = aJ(A;i,0) _1 = — ik\ which corresponds exactly to 
the standard action with nearest neighbor couplings. This is maximally local. For 
general 2-dimensional configurations the fixed point action with c = 9 not only 
contains nearest neighbor interactions. However, as shown in fig. 2 the corresponding 
function p\{z) in coordinate space is still very local. Some values of p\{z) are given in 
table 1. Being a quantity with a tilde pX(z\, z 2 ) refers to the shifted fermions located 
at the sites of the lattice with spacing 1/2. It describes the coupling between the 
pseudoflavors 1 and 2, which are separated by a half-odd integer multiple of the 
lattice spacing in the 1-direction and by an integer multiple of the lattice spacing 
in the 2-direction. Hence z\ is a half-odd integer while z 2 is an integer. This 
also explains the boundary conditions of the corresponding quantity p\{k) (and 
similarly a\(k)) in momentum space. In the 1-direction both are antiperiodic over 
the Brillouin zone due to the corresponding half-odd integer differences in coordinate 
space, and they are periodic in the 2-direction because there the differences are 
integers. 

3 Fixed point actions for the Gross-Neveu model 
in the small field approximation 

In this section we turn to the Gross-Neveu model [8], i.e. we switch on a 4- Fermi 
interaction which is linearized by a Yukawa coupling to a real valued auxiliary scalar 
field. The interaction breaks the U(l) e <E> U(l) chiral symmetry of the free theory 
explicitly down to [/(l) e=0 ®Z(2). Still, the Z(2) symmetry is sufficient to protect the 
fermions against perturbative radiative mass corrections. However, the interaction 



11 







1 



2 3 
z 



4 



Figure 2: Exponential decay of the coupling \f>l(\z\)\ as a function of \z\ for the 
optimally local fixed point action with c = 9. The triangles correspond to difference 
vectors parallel to the lattice axes. The straight line is a fit to the triangles resulting 
in 1^(^,0)1 oc exp(— 3.46|zi|). 

is strong enough to break the Z(2) symmetry spontaneously which results in a 
nonperturbatively generated fermion mass. The corresponding order parameter is 
the chiral condensate. 

As it was shown in ref. [10] the auxiliary scalar held $ should be placed at the 
centers of the plaquettes of the lattice with spacing 1 /2. This geometry — consistent 
with the so-called hypercubic Yukawa coupling — is illustrated in fig. 3. We now 
parametrize the action as 

^S[x, X , $] = Xxp{x - y)Xy + E \®l + E Xx°{x -z,y- z) Xy $ z }. (3.1) 

Here a(v,w) is a 4 X 4 matrix in pseudoflavor space that depends on the differ- 
ence vectors v = x — z and w = y — z between the fermionic and auxiliary scalar 
variables. Again the symmetries of the problem impose constraints. The remnant 
chiral symmetry implies cri 2 (v } w) = cr 2 i(v } w) = 0, cri 3 (v } w) = cr 3 i(v } w) = 0, 
cr 2 4(v } w) = o-4 2 (v } w) = 0, <734 (u,u;) = o-4 3 (v } w) = 0, charge conjugation leads to 
aij(w,v) = aji(v,w), and the shift symmetries impose further constraints that we 
don't write down explicitly. 

To determine the fixed point action for the interacting theory one must also 
define a renormalization group transformation for the auxiliary scalar held. Since 
the scalar held does not carry a pseudoflavor index one could form ordinary 3x3 
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Z\ 


Z2 


P\{zi, z 2 ) 


0.5 


0.0 


0.8656403 


1.5 


0.0 


0.0161216 


2.5 


0.0 


0.0005065 


3.5 


0.0 


0.0000174 


0.5 


1.0 


0.0649779 


1.5 


1.0 


-0.0072205 


2.5 


1.0 


-0.0003117 


0.5 


2.0 


0.0021230 


1.5 


2.0 


-0.0007996 


0.5 


3.0 


0.0000788 



Table 1: Some values of the function p\{z) for the optimally local fixed point action 
with c = 9. 

blocks of neighboring variables. Here we prefer to work with scalar field blocks that 
are of the same kind as the fermionic ones. This is illustrated in fig. 3. One may even 
argue that a scalar pseudoflavor is defined by the surrounding fermionic variables 
on the corresponding plaquette. Blocking schemes that respect this structure are 
limited to blocking factors 5, 9, 13, etc. We have checked explicitly that the resulting 
fixed point actions are identical with the ones we derive here. 

First, we determine the fixed point action for the free auxiliary field. The renor- 
malization group step is then defined as 

exp(-4<7'[$']) = /^exp(-i^i^ )exp( _^^ ( ^,_A^^)2). (3.2) 

z z 1 z£z' 

Finding the fixed point is much simpler here than it was in the fermionic case. In 
fact, one need not even go to momentum space. We simply quote the result 

S '^ = ^ a -¥T^^- < 3 - 3 > 

Thus the fixed point condition reads 

The auxiliary field remains auxiliary after the renormalization group step, in par- 
ticular, no kinetic term is generated. Thus, the dimension of $ is o?$ = 1 and not 
zero as for a genuine dynamical scalar field in two dimensions. Hence, we expect 
fS = 3 d * = 3. Then eq.(3.4) implies a = oo, i.e. we are restricted to a ^-function 
renormalization group transformation for the auxiliary field. We will see later that 
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Figure 3: Geometry of the fermionic and auxiliary scalar variables in the Gross- 
Neveu model with hypercubic Yukawa coupling. The crosses represent the auxiliary 
scalar field. A 3 X 3 block of auxiliary field variables is also shown. The nine crosses 
that contribute to the block variable are marked with a small circle. 

this is indeed the case. For the moment we leave a and f3 undetermined. When we 
come to the interacting theory we need the relation between the classical solution 
$ c and the blocked field $' in momentum space which is given by 

$ c (*) = ^-^mWik). (3.5) 

It should be noted that the fixed point action for the free auxiliary field does not 
correspond to a critical theory. In fact, the correlation length of the auxiliary field 
vanishes instead of being infinite. 

In the next step we determine the fixed point action for the interacting theory 
in the small field approximation, i.e. we assume that $ is small and we restrict 
ourselves to respectively Up to now we have determined the blocked 

action as 

s'\x', x', = s'\x', x'] + SW] + om. (3.6) 

Since we are at G — > the value of S'[x', x'] was determined by inserting the solutions 
X c , Xc of classical equations of motion (2.19,2.20) into the exponent of eq.(2.12). We 
write symbolically 

S'[x' } x'] = E[ Xc ,Xc]. (3.7) 

Similarly S 7 ^'] = _E[$ C ]. The exponent we are dealing with now contains the 
Yukawa interaction as an additional term of 0($) 

E[ X , X, $] = E[ X , X] + E[$] + Y[ X , X, (3.8) 

This term changes the classical equations of motion and hence their solutions at 
0($') to Xc + SXc + 0($ /2 ), Xc + SXc + 0($ /2 ), and $ c + <5$ c + 0($ /2 ), where S Xc , 
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Sxc and <5$ c are of To determine the blocked action to one must insert 

the new solutions into the new exponent 

s"[x',x',$'] = % + %,Xc + %,$c + ^ c ] + o($' 2 ) 

= E[xc + 6 Xc , X c + S X c] + E[<$> c + 5$ c ] + y [ Xc , Xc , $ c ] + 0($' 2 ) 

= E[ Xc , Xc] + £[$c] + y [Xc, Xc $c] + 0($' 2 ) 

= S'[x\x'] + S'm + Y[x cl Xc^c] + 0(^). (3.9) 

The essential observation is that Xc, Xc an d ^c are solutions of classical equations of 
motion, i.e. extrema of the corresponding exponents. This implies E[x c + 8x c , Xc + 
Sxc] = E[ Xc ,Xc] + 0($' 2 ) and E[$ c + «5$ c ] = £[$ c ] + 0($' 2 ). Consequently, to 
determine the blocked action to (9($ / ) at G — > it is sufficient to insert the old 
solutions of the equations of motion into the Yukawa term. For this purpose we 
again go to momentum space. Then the Yukawa term takes the form 

Y[x, x, $] = 7A7 /_ d 2 pd 2 qx{-p)v{p,q)x(-q)${p + q)- (3.10) 

{iir) jb 2 



The integration of both momenta p and q extends over a large Brillouin zone B = 
] — 27r, 2tt] 2 because the auxiliary scalar field lives on a lattice with spacing 1/2. To 
obtain the Yukawa term after one renormalization group step we simply insert the 
old solutions of the classical equations of motion from eqs. (2. 19,2.20,3.5) into this 
expression and we find 

Y'[x',x',&] = Y[xc,Xc$c\ 

,2 32 «/3 1 



32 + a/3 2 ( 47r )4 J m d 2 pd 2 qx(-p)a(3p)- 1 U(p)D(p)a(p) 
X <r(p, q)a{-q)D{q)Il{q)a , {-3q)- 1 X , {-q)Il{p + q)&(p + q) 

■pj I d 2 pd 2 qx'(-py(3p,3q)x'(-qW(p + q)- (3-11) 

47T/ JB 12 

After the renormalization group step the integration extends over the smaller Bril- 
louin zone B' =] — 27r/3, 27r/3] 2 . One reads off the recursion relation 

U ' H > 3 2 + a/3 2 3 4 f^ UJ y 3 J y 3 ; 

p + Airl q + i-Kin q + Atttti q + 47rm , 
x CT ( 3 , 3 )«( 3 )D{ )«(-?) 

„ , p + 47r/ , „ , o 4- 47rm , „ , p 4- 47r/ 4- q + 47rm , 

x n(^— )n(^ — )n(£± — ). (3.12) 

Introducing a(p,q) = D(—p) 1 / 2 a(p)(7(p,q)a(—q)D(—q) 1 / 2 , eq.(3.12) turns into the 
recursion relation 

2 a/3 1 p + 47t/ 5 + 47rm 
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ri(-i) 



,„ +m „ 2sin(p„/2) 2sin( g „/2) 

6sin((p i/ + 47r/ J/ )/6) 6sin(( gi/ + 47rm J/ )/6) 
2sin((p J/ + ^)/2) 



6sin((p i/ + 47r/j, + (jv + 47rm i/ )/6) 
Iterating this we obtain after n renormalization group steps 



(3.13) 



a/3 



1 



V 3 2 + a/3 2 3 2 



E 



,p + 47r/ ^ + 47rm, 



2sin(^/2) 



2sin(^/2) 



x 



3 n 2 sin((^ + \tt\ v ) /3 n 2) 3 n 2 sin(( gl/ + 4^,,) /3"2) 
2sin((^ + ^)/2) 



3 n 2 sin((^ + \tt\ v + q„ + 4 7 rm i ,)/>2) 
The standard action is characterized by 



(3.14) 



a{p,q) 



( *o(p,q) 

a (p, q) a 3 (p, q) 

-a 3 (p,q) cr (p,q) 

V Hp, ?) 



-°3(p,q) \ 


o 



(3.15) 



with 



P^u 
p 2 q 2 



n 



tij.pPij.qp 



cos 



p g 4 



,£*/ + qv 



4 y ' 

+ qv 



(3.16) 



Here e^p is the antisymmetric tensor in two dimensions. In the limit of infinitely 
many renormalization group steps only the behavior of <r(p, q) close to its poles 
matters. The poles that contribute to the fixed point action are at p = q = 2iri with 
ifj, £ {0, 1}. After infinitely many steps one obtains 



K(p,q) 



a/3 



Km U % 2 , R2 



EE 



(^ + 471-/^ + 2^i jl ){q jl + Airm^ + 2^) 



n(-i) 



^ m t . (p + 47r/ + 27Tz) 2 ( g + 47rm + 27rz) 5 
;, +m , +8 , 2sin(p„/2) 2sin( g „/2) 

p v + 47t/j, + 27TZJ, q v + 4^77^ + 27TZJ, 

2sin((^ + q„)/2) 



p v + 4^/^ + q v + 4^771^ + 47TZj, 
A nontrivial fixed point is reached only if 

7.2 a P 



3 2 + a/3 2 



1. 



(3.17) 



(3.18) 



Using b 2 = 3 together with the fixed point condition eq.(3.4) for the free auxiliary 
scalar field this implies a = oo and (3 = 3. Hence, one is limited to a ^-function 
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renormalization group transformation for the auxiliary field. The canonical dimen- 
sion of the auxiliary field is o?$ = 1 consistent with fi = 3 d * = 3 as discussed earlier. 
Only for this value the Yukawa coupling in two dimensions is a marginal operator. 
In particular, for a genuine dynamical scalar field (which would have o?$ = in two 
dimensions) the Yukawa coupling would be irrelevant. In fact, the Yukawa coupling 
to the auxiliary scalar field is weakly relevant. This effect is, however, invisible in 
the leading order of the weak field approximation. All this is expected because just 
in two dimensions the Gross-Neveu model is asymptotically free. At the fixed point 
we finally obtain 

s (P» + 4?r ^ + 27rz At )(5 At + Airm^ + 2iri ll ) 



2 • (p + 47r/ + 2Tri) 2 (q + 47rm + 2tti) 



jt, 1 y y+nty +,- y 2sin(p„/2) 2sin(g„/2) 
2sm((p l/ + q l/ )/2) 



p v + AttIv + q v + AttuIv + \ki v 

(Pf, + Airlf, + 2iri ll )(q p + 47rm p + 2tti p ) 



v2 . (p + 4tt/ + 2ttz) 2 (<7 + 4vrm + 2ttz) 2 
tt/ ^L+my 2sin(p„/2) 2sin(g„/2) 

2sm((p l/ + q l/ )/2) 



p v + 47t/j, + q v + 4^77^ + 47TZj, 



(3.19) 



These functions are 47r-antiperiodic in both momenta over the large Brillouin zone 
B. Being quantities with a tilde, they refer to the shifted fermions living on the sites 
of the lattice with spacing 1/2. Then the fermionic and auxiliary scalar degrees of 
freedom are separated by multiples of 1/4 of the lattice spacing. In momentum 
space this would imply ±z-periodic boundary conditions over the ordinary Brillouin 
zone 73, and it implies antiperiodic boundary conditions over the larger Brillouin 
zone B. 



4 Fixed point actions for the Gross-Neveu model 
at large N 

In the 0(3) model Hasenfratz and Niedermayer have determined the full (large field) 
fixed point action using a numerical minimization procedure on a multigrid [6]. Still, 
in their calculation the small field approximation was essential because it allowed 
them to optimize the fixed point action's locality analytically. In general it will be 
impossible to go beyond the small field approximation using analytic methods. The 
only exceptions we are aware of are free field theories and theories in the large N 
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limit. Here we consider the Gross-Neveu model, but the method should work just 
as well e.g. for 0(N) and CP(N) models. For N = oo the fixed point of the O(N) 
model has been studied in ref. [15]. For a large number N of flavors the dynamics of 
the Gross-Neveu model greatly simplifies, because then the auxiliary scalar field is 
dominated by a constant mean field $ , whose value is determined by the minimum 
of its effective potential. In fact, in the large N limit only the constant zero mode 
can have a large value while the amplitudes of the nonzero modes are suppressed by 
powers of 1/N and are thus small. Hence, the small field approximation is justified 
for the nonzero modes, and only the zero mode needs special treatment. First we 
concentrate entirely on the zero mode assuming that the auxiliary scalar field is a 
constant $ - Later we include the nonzero modes using the small field approximation 
to leading order. Since a constant field remains constant after ^-function blocking, 
the renormalization group step for the auxiliary field reduces to 

& = /3$ = 3$ . (4.1) 

From now on we consider N/2 sets of staggered fermions (with N even) correspond- 
ing to N flavors of Dirac fermions in the continuum limit. We suppress the additional 
flavor index. In the following expressions summation over all flavors is understood. 
In momentum space the constant auxiliary field is $(&) = (47r) 2 $ ^(^) and the 
Yukawa term takes the form 

y[x,X,$o] = 77^ [_d 2 kx(-k)a(k,-k)x(k)% 
(4^)^ Jb 

1 d 2 k x (-k)X(k)l x (k). (4.2) 



(2vr) 2 

Here we have used that x{ — k) and x(k) are 27r-periodic and we have defined 

X(k)l = ^£>(& + 2™,-A;-2™)$o, (4.3) 

i 

where the sum extends over integer vectors with components £ {0,1}. For the 
standard action one finds X(k) = $ - It takes the form 

S\x,x,* ] = i^-r 2 I d 2 kx(-k)[p(k) + x(k)l}x(k) + ^lv 

{Air) -Jb 2 
= j^j B d 2 kx(-k)[p%k) + %l}x(k) + UlV, (4.4) 

where V is the space-time volume. The Yukawa term acts as a mass term of the 
fermion field. To find the fixed point action we proceed as in section 2. The calcula- 
tion is exactly the same except that we now define a(k) + f3(k)1 = [p(k) + X(k)l.]~ 1 . 
In particular, the recursion relation eq.(2.23) remains unchanged except that it now 
holds for a(k) + f3(k)1. (We do not introduce f3(k) because it would be identical 
with f3(k).) Now the initial matrix elements characterizing the standard action are 
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given by a^k) ~ ik^j ' {k 2 + $q) and f3(k) ~ $ /(^ 2 + $0) f° r sma ll momenta. Thus, 
at the fixed point one finds 

Next we compare this result with the weak field approximation of the previous 
section. To (9($ ), a *^{k) is identical with the one of the free theory in eq.(2.26). 
We define \*{k)l = &*(k)X*(k)l&*(k) and we find 

a-,,, = - nk) + o«> = - e j^tr n (xrlf ) 2 + °«>- <«» 

On the other hand, analogous to eq.(4.3) one expects 

\*{k)l = ^J2 ^*( k + 2 ™> ~* " 2 ™)$o. (4-7) 
From eq.(3.19) one obtains ^(k } —k) = and 

a ° [ ' ' (k + iirl + 2in) 2 (-k + iirm + 2tti) 2 

x \\{ piv+mv+iy 2sin (^/ 2 ) 2sin(-^/2) 

V ^ + 47T/, + 27TZ, -k v + 4^, + 27TZ, 

1 yrf 2shl(^/2) V 

^ 2 Y {k + 4vr/ + 2vrz) 2 \ [ {K + Airl + 2« J 

= v 1 n f 2sin( ^ /2) V a*) 

^ (k + 2vr/)2 A/ ^ + 2^ y 1 ' ; 

such that indeed 

ly>*(£ + 27rz,-£-27rz) = - V —J—rrf^^Vl. (4.9) 

In the next step we include the fluctuations of the nonzero modes. In the large N 
limit they are small and can be treated by the small field approximation. Later we 
will be interested only in the fermion mass and in the chiral condensate. To leading 
order in 1/N these quantities don't pick up contributions from the nonzero modes. 
However, when one wants to compute e.g. the mass of the fermion-antifermion 
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bound state (which we don't do in this paper), the dynamics of the nonzero modes 
enters at the lowest nontrivial order in 1/N. We parametrize the Yukawa coupling as 
in eq.(3.10). However, now does not contain the zero mode $ that was taken 

into account separately. The calculation of the fixed point action is very similar to 
the one in section 3. The only difference is that now a(p,q) = D(—p) 1 ^ 2 [a(p) + 
f3(p)\\a(p,q)[a(—q) + f3(— q)l]D(— q) 1 ^ 2 which has the structure 



a{p,q) 



( °o(p,q) 
^(p,q) 



°i(p,q) °2(p,q) 

&o(p, q) cn{p, q) 

-a 3 (p,q) v (p,q) 

-cr 2 (p,?) cri{p,q) 



-°3(p,q) \ 
-^(p,q) 
&i(p,q) 
&o(p, q) 



(4-10) 



The recursion is now started with 

&o(p, q) 
°n(p,q) 
&3(p, q) 



p„q„ + $o 



(p 2 + $8)(? 2 + $8; 

ih - q^o 



(p 2 + $g)(? 2 + $t 

tij.pPij.qp 



n 

V 



cos 



COS 



(p 2 + $g)(? 2 + K 



n 



COS 



,p v + q v 

4 • 

,p v + qv 
4 

,p v + q v 



(4-11) 



At the fixed point this leads to 



(J, 



o(p,q) 



E E 



fa + 4?^ + 2 7 rz At )(g At + Awm^ + 2?r^) + $ 2 
[(p + 4tt/ + 2ttz) 2 + $§][(? + 4vrm + 2ttz) 2 + $ 2 



n(-i) 



2sinfa/2) 



2sin(^/2) 



+ 4^/^ + 2iri u ^ + 4^771^ + 2iri l 
2sm((p l/ + q l/ )/2) 



K{p,q) 



p v + 47r/j, + (JV + 4^ 777^ + 47T^ 

fa + Airlf, - 



E E 



^ -4™,)$ (-l) ! '" 



[(p + 47r/ + 27TZ) 2 + + 4vrm + 27rZ )2 + $2] 

2sinfa/2) 2sin(^/2) 



p v + 47t/j, + 27TZj, ^ + 4^ 777^ + 27TZj, 

2sin(fa + ^)/2) 



+ 47r/j, + (JV + 4^ 777^ + 4:7V i L 



°t{p,q) 



E E 



l.ra 



tijpjPp. + 4^/^ + 27ri ll )(q p + 47rm p + 2txi p ) 
[(p + 4tt/ + 2ttz) 2 + + 47rm + 2ttz) 2 + $g] 

2sinfa/2) 2sin(^/2) 



p v + 4^/^ + 2%i v q v + 4^777,^ + 2%i v 
2sin(fa + ^)/2) 



+ 47r/j, + (JV + 4^ 777^ + 4^7^ 



(4-12) 
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Finally, the fixed point action is given by 



+ 77^7 /- d 2 pd 2 qx(-p)a*(p,q)x(-qMp + ?)}. (4.13) 
(47rJ^ jb 2 

Only the first two terms contribute at N = oo. The other terms represent the 
leading order of the small field approximation for the nonzero modes. They are 
important only at higher orders in 1/N. It should be noted that the small field 
expansion can be extended systematically to higher orders. The above expression 
for the fixed point action has been derived in the limit G — > 0. We like to point out 
that it is the analog of the perfect classical action of Hasenfratz and Niedermayer 
that they determined numerically in the 0(3) model. Finding the fixed point action 
is a saddle point problem. In fact, we have found the fixed point just by solving 
classical equations of motion. Although we used a path integral formulation, so 
far we have only looked for saddle points. In this sense the fixed point action is a 
classical action. Of course, the path integrals one encounters in the large N limit are 
all saddle point problems, and solving the classical problem is equivalent to doing 
the full integral. This is the reason why the fixed point action is even a perfect 
quantum action, i.e. there are no cut-off effects even if we leave the critical surface. 
This is the case only in the large N limit. For finite N one would expect some cut-off 
effects at finite correlation lengths even with the fixed point action. Still, the fact 
that the perfect classical action is a perfect quantum action in the large N limit 
may explain why it is practically perfect in the 0(3) model studied by Hasenfratz 
and Niedermayer. 

Next we construct a perfect operator for the chiral condensate. The standard 
operator is given by 

Ex«Xx = j^y 2 J B d 2 kx(-k) X (k). (4.14) 

The chiral condensate is a relevant operator — it gets amplified under renormal- 
ization group transformations and hence it is driven away from the fixed point. To 
determine the perfect operator we introduce a small perturbation around the fixed 
point 

SAx, X, $] = s*[ x , X, $] + jX[ X , X, (4-15) 

The perfect operator is an eigenfunctional that reproduces itself under 

renormalization, i.e. 

^x',^] = s*' [x\x j x*'[x\x = s*[x\x'^']+nx*[x\x'^']+0(j 2 ). 

(4.16) 
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For relevant operators 7 > 1. In two dimensions the chiral condensate has dimension 
dy? x = 1. Hence, the corresponding eigenvalue is 7 = 3 dxx = 3. Later we will be 
interested in the vacuum value of the perfect operator. To leading order in 1/N 
this quantity gets contributions only from the zero mode of the scalar auxiliary 
field. Hence we may assume that the scalar field is a constant $ - Still, if one 
wants to compute e.g. the perfect correlation function of \X a ^ so the nonzero modes 
contribute at the lowest nontrivial order in 1/N. It is straightforward (but somewhat 
tedious) to derive the nonzero mode contribution in the small field approximation. 
Here we restrict ourselves to the zero mode contribution to the chiral condensate 
which we parametrize as 

X[X, X, %] = j^y j B d 2 k x(-k)Hk) + u(k)l] X (k). (4.17) 

Here fi(k) is a matrix of the same structure as p(k). In general both fi(k) and 
v(k) are functions of $ - The standard operator is characterized by fi(k) = and 
v(k) = 1. To determine the perfect operator we proceed analogous to the perfect 
action. First we define [a(k) 3 + ^(k)^!]' 1 = p(k) + X(k)l + j[p(k) + v(k)l] as 
well as a(k)j = D(-k) 1 l 2 a(k) ] D(k) 1 l' 2 and fi(k) = D(-k) 1 / 2 /i(k)D(k)^ 2 . Then the 
recursion relation eq.(2.23) holds for ct(k)j + f3(k)j1, except that now 

[«'(% + P'Wjl]- 1 = p'(k) + \'(k)l + m'{k) + u'(k)l] + 0(j 2 ). (4.18) 

Iterating eq.(2.23) one finds after infinitely many steps 

&*M = &;(k) + J Hm 3 n d* &l(k) + O(j 2 ), 

ff{k) 3 = /3*(k)+j]im o rd 9o l3*(k) + 0(j 2 ). (4.19) 

This confirms 7 = 3 and it allows to identify the perfect operator that is character- 
ized by 



= a 



$0 



&t(-k)&t(k) + p*(-k)p*(k). 



Note that a*(k) is odd in k^ and even in the other components k v ,v ^ //, while 
f3*(k) is even in all components of the momentum. 

How local is the perfect operator for the chiral condensate? Here we answer this 
question only in the small field approximation omitting O($o) contributions, and 
we consider again effectively 1-dimensional configurations. From eq.(4.5) we obtain 
f3*(k) = O($ ), d$ al(k) = O($ ) and 

iteZ ( 1 + 2irJ 'V \ h + 2tt/i J k\ § 
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Using eq.(4.20) together with eq.(2.27) for c = 9 one finds p*^(k) = 0($o) an d 

v*(k 1 ,0) = l-h 2 1 + O($ 2 ). (4.22) 



In coordinate space this implies for the perfect chiral condensate in the small field 
approximation 

ixx)x = \,XxXx + -^{xxx x+ i + x x+ iXx + XxX x -i + x x -iX*), ( 4 -23) 

which is extremely local. For general 2-dimensional configurations the perfect chiral 
condensate not only contains nearest neighbor contributions, but still it is very local. 



5 Solving the Gross-Neveu model at large N 

In this section we solve the Gross-Neveu model at large N using the full and small 
field fixed point actions as well as the standard action. For an explanation of the 
large N technique we refer to Coleman's book [16]. We compute the energy spectrum 
of fermionic 1-particle states and the chiral condensate. This shows that the full fixed 
point action is a perfect quantum action. Up to now we were always on the critical 
surface, i.e. at infinite correlation length. Now we will use the perfect classical (fixed 
point) action away from the critical surface in a region where quantum fluctuations 
are present. 

To solve the model we first compute the effective potential V^//($o) f° r constant 
auxiliary scalar fields $ 

exp(-U eJJ ($ )U) = J V X V X ex ? (-^S[ X ,xM) = detA^o)^ 2 exp(-^U), 

(5-1) 

where V is the space-time volume, Ai(<&o) is the fermion matrix of a single set of 
staggered fermions in the constant background field $ , and the power N/2 accounts 
for the various flavors. Up to now S[x, X? ^o] is still general. We parametrize it as 
before 

S[x,X^o] = J7 ^- 2 j d 2 kx(-k)[p(k) + X(k)l]x(k) + U 2 V } (5.2) 
{Air) JB 2 

where p(k) and X(k) are functions of $ - One obtains 

Veff(%) = t^K ~ f j^y 2 j B d 2 klogdet^[p(k) + X(k)l]. (5.3) 
The minimum of the effective potential is given by 

cKV eff (%) = -- -j^y 2 J b d k det[m + mi] = 0. (5.4) 
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Here we have used det [/?(&) + A(&)1] = det [/?(&) + A(&)1]. Eq.(5.4) is an implicit 
equation for $ , known as the gap equation. It determines the relation between the 
bare coupling G and the vacuum value of the auxiliary field $ - For the standard 
action p^ik) = —ik^ and X(k) = $ - Then the gap equation reduces to 

1 1 d 2 k^ . (5.5) 



g (2tt) 2 Jb k 2 + $g 

We recall that the large N limit is taken such that the coupling g = GN is kept 
fixed. The gap equation for the fixed point action is given by 

$o_ 1 f £k 2[&t(-k)d 9o at(k) + P*(-k)d* /3*(k)] 

~ -te. I 7 _ \ ~ ale / 7 . \ I Q-te. I 7_ \ /Oilt / 7_ \ V / 



g (2vr) 2 Jb a*(-A;)a*(A;) + 0*(-k)f3*(k) 

Next we investigate the spectrum of the fermionic 1-particle states. For purely 
technical reasons related to the proper normalization of momentum eigenstates the 
system is put in a finite spatial volume of size L with periodic boundary conditions. 
We consider the correlation function of a fermionic operator with definite spatial 
momentum k\ = 27rni/ ' L £] — ~, ~], // 1 £ Z 

1 f n 

x(ki)x 2 = — dk2x(h,h)exp(ik 2 X2) (5.7) 

Z7T J-tt 

that creates a fermion at euclidean time x 2 . The correlation function is given by 

L r N 

{x(-ki)ox(ki) X2 ) = TT / ^ 2 — TrA^($ )" 1 exp(iA; 2a;2 ). (5.8) 

Z7T J-tt Z 

Here the trace is over the pseudoflavor index of a single set of staggered fermions 
only. The inverse fermion matrix is given by ^(^o) -1 = —G[a(k) + 0(k)1]. For 
the standard action we have 0(k) = $ /(^ 2 + ^o) sucn that 

(x(-h)ox(hU) = f dk 2 2 f ° e M ik ^ = Cik^eM-Eih^). (5.9) 

27T J — TT k\ + k\ + $o 

As expected the correlation function decays exponentially, and the energy of the 
1-fermion state is given by the pole of the integrand as E(k\) = —ik 2 with 

4sinh 2 (£(£4)/2) = -4 sin 2 (& 2 /2) = -k\ = k\ + ( 5-10 ) 

The fermion mass mj is the energy of the k\ = state, i.e. 

2sinh(m / /2) = $„• (5.11) 

For small mj and k\ the momentum dispersion relation turns into the one of the 
continuum theory. However, there are 0(a 2 ) corrections to it. It is interesting to 
note that all cut-off effects are 0(a 2 ) } while one would naively expect them to be of 
0(a) in a fermionic model. Since we have not introduced an explicit chiral symmetry 
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breaking fermion mass term, the action is Z(2) chirally invariant. This symmetry 
prevents the generation of irrelevant operators of dimension 3, which would induce 
0(a) scaling violations. For the same reason the standard action is in this case 
equivalent to an 0(a) Symanzik improved action. Hence, it is of the same quality 
as the Sheikholeslami-Wohlert action for QCD. For the small field fixed point action 
one finds 



1 - Jl - 2$2/ 3 
2 sinli(m / /2) = , (5.12) 

and again there are 0(a 2 ) corrections to the continuum dispersion relation. The 
small field fixed point action is definitely not a perfect action. Finally we consider 
the full fixed point action. Then using eq.(4.5) the eq.(5.8) turns into 

L r ji 2#$ exp(z& 2 £ 2 ) tt f 2 sm (W 2 )\ 2 



^ (k + 2tt/) 2 + $g 11 y kiy + 27r l J 
L f°° ^ 2g$ exp(ik 2 x 2 ) (2 sin(&i/2) \ 2 (2 sin(k 2 /2) \ 2 



(x(-h)ox(hu) = -- j_ v dk 2 e (k+27rlY+ ^ n n 

/oo 
dk 
-00 2 ^ (h + 2Trh) 2 + k\ + $ 2 \k 1 + 2tt/i J \ k 2 J 

= C ( k ^ + 2?r/ i) exp(- J E(A; 1 + 2 7 r/ 1 )x 2 ). (5.13) 

The sum over l 2 and the integration of k 2 over the Brillouin zone have combined to 
an integration over the momentum space of the continuum theory. The remaining 
summation over li leads to infinitely many poles of the integrand, and hence to 
infinitely many states that contribute to the correlation function. Their energies are 
given by E(k\ + 2tt!i) = —ik 2 with 

E 2 (k 1 + 2wh) = -k 2 = (h + 2tt1 1 ) 2 + (5.14) 

This is the exact continuum dispersion relation for a particle of mass m/ = $ 
with momentum k\ + 2 r xl\. Remarkably, this momentum is unlimited although the 
lattice momentum k\ is restricted to the Brillouin zone ] — 7r,7r]. This means that 
the full fixed point action has exactly the same fermionic 1-particle spectrum as 
the continuum theory. There are no cut-off effects — the action is perfect. It is 
interesting to see how the perfect action can account for all momentum states of the 
continuum theory. Since it is a lattice action it is invariant only against translations 
by multiples of the lattice spacing. Hence, like for any other lattice action, using 
the discrete translation symmetry one can only select states which have a definite 
lattice momentum k\ G] — 7r, tt]. Still, for the fixed point action the other momentum 
states k\ + 2tt!i of the continuum theory show up, and they all contribute to the 
correlation function with the lattice momentum k\ that was selected by symmetry. 
This is possible only because the fixed point action is extended in euclidean time. In 
principle, it extends over all time slices, but it falls off extremely fast, and it is hence 
practically local. Still, if one wants to catch even the highest states of the continuum 
theory, one must take the coupling of all time slices into account. In practice, e.g. 
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in a Monte-Carlo simulation, one would not like to keep the exponentially small 
tails of the fixed point action. When one cuts them off one does some harm to the 
continuum spectrum. However, this only affects the highest states, while the energies 
of low lying states are still practically perfect. The same remarks apply to rotation 
invariance. As a lattice action the perfect action is not spherically symmetric. Its 
symmetry is restricted to the discrete cubic rotation group. Still, the spectrum does 
not notice this — it is the exact one of the continuum theory. The lattice symmetry 
only allows us to construct operators that project on a given representation of the 
cubic rotation group. Still, all states of the continuum theory have some projection 
on these operators. Remarkably, when the fixed point action is used, their energy is 
perfect. Again, this is possible only because the action is extended in space. Cutting 
off its exponentially suppressed tail will do some harm, but it will mostly affect the 
high angular momentum states that one does not focus on in numerical studies. 

Now let us consider the chiral condensate. Using the standard action together 
with the standard operator for the chiral condensate one finds 

(xx) = j^y 2 I d 2 k ^TrM(^)- 1 = I d*k J^- = -$ . (5.15) 

{Ztt) 2 Jb 2 (2tt) 2 Jb k 2 + $g 

Here we have used the gap equation (5.5). The ratio {xx)/ m f should scale, i.e. it 
should be cut-off independent in the continuum limit. Indeed 

{xx)/ m f = ~~ 2 sinh(mj/2)/mj, (5.16) 

which approaches the constant —I in the continuum limit £ = l/mj — > oo. However, 
there are 0(a 2 ) cut-off effects as long as the correlation length £ is finite. Using the 
perfect action with the perfect operator for the chiral condensate one obtains 



(XX) = -TTT^i / <PklTr{\pL*(k) + v*(k)l][a*(k) + P*(k)l]} 



(2tt) 2 Jb 2 

i / * h 2g[&:(-k)d* &:(k) + p*(-k)d* p*(k)] _ f ^ 



(2tt) 2 Jb &l(-k)&l(k) + f3*(-k)/3*(k) 

Here we have used eqs. (4.20,5.6). Now {xx)/ m f = — 1 even for arbitrarily small £. 
Thus, cut-off effects are completely eliminated. It should be noted that this does 
not happen automatically just by using the perfect action. Also the operator for the 
chiral condensate must be perfect. With the standard operator XxXx there would 
be cut-off effects even when the perfect action is used. 

Finally, we investigate the asymptotic scaling behavior. For the standard action 
this has already been done in ref.[ll]. Fig. 4 depicts rrif(g) for the perfect action and 
for the standard action. For the perfect action asymptotic scaling sets in earlier than 
for the standard action. In fact, for the perfect action the deviations from asymptotic 
scaling are only about 5 % for correlation lengths as small as 1. Note that this 
corresponds to 2 lattice units on the fine lattice (with spacing 1/2). Still, asymptotic 
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1/g 

Figure 4: Scaling behavior rrif(g) of the perfect action (solid line) and of the stan- 
dard action (dashed line). Asymptotic scaling with the perturbative f3 -function cor- 
responds to the straight dotted lines. 

scaling is not perfect for the perfect action. Fig. 5 compares the /3-function /3(g) = 
—mfd mj g of the perfect and the standard action with the perturbative one obtained 
from eq.(1.5) 

/%) = -9 2 /*. (5-18) 

Again, one sees that asymptotic scaling sets in earlier with the perfect action, but 
it is not perfect. Asymptotic scaling is not what one should ask for. What really 
matters is that the physics (dimensionless ratios of quantities in lattice units) are 
cut-off independent. This is the case for the perfect action. The fact that it does not 
show perfect asymptotic scaling is unphysical because it involves the bare coupling 
defined at the cut-off scale. Scaling, on the other hand, is a physical issue, and 
scaling is perfect for the fixed point action. 

It may still seem mysterious why the fixed point action — which is a perfect 
classical action — is even perfect on the quantum level, i.e. why there are no cut- 
off effects even at arbitrarily small correlation length. To shed some light on this 
question we investigate the renormalized trajectory of theories with perfect quantum 
actions in the large N limit. We want to show that the classical perfect action is 
in fact located on the trajectory. Let us say we are aiming at the perfect quantum 
action for the theory at some finite correlation length £ = 1/mj. We construct 
this action by performing n exact renormalization group steps starting at a point 
arbitrarily close to (but not on) the critical surface, i.e. we start at a correlation 
length 3 n ^. We are interested in the limit n — > oo with mj = l/£ fixed. Of course, 



27 




0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 

Figure 5: The /3-function (more precisely /3(g) / g 2 versus l/g) for the perfect action 
(solid line) and for the standard action (dashed line). Asymptotic scaling with the 
perturbative /3-function corresponds to the dotted line. 

now we are at g ^ and for finite N performing a renormalization group step is 
no longer a saddle point problem. Instead one must perform the full path integral. 
At large A, however, we may again assume that only the zero mode $ of the 
auxiliary scalar held has a large value, while the nonzero modes can be treated 
in the small held approximation. Then again all path integrals are saddle point 
problems, and their result is identical with the classical value. Thus, an infinite 
number of renormalization group steps leading us down the renormalized trajectory 
results again in the perfect classical action, except that now g ^ 0. Clearly, this 
is the case only in the large A limit. At finite A one expects the perfect classical 
action to have some cut-off effects at small correlation lengths. 

6 Conclusions 

The main purpose of the present paper was to investigate fermionic fixed point 
actions in a simple setting, and to understand better the success of Hasenfratz's 
and Niedermayer's perfect action approach in an analytic manner. This confined us 
to the large A limit, because only there we could construct the fixed point action 
analytically. It turned out that the perfect classical (fixed point) action is in fact 
also a perfect quantum action. Cut-off effects are completely eliminated even at 
arbitrarily small correlation lengths. Certainly, this is the case only in the large A 
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limit. It is related to the fact that then all path integrals are saddle point problems 
and classical methods already give the full answer. Still, we believe that this may 
explain why the classical fixed point action worked so well in the 0(3) model. The 
fact that the perfect classical action is perfect on the quantum level at large N 
may imply that it is still very good for N = 3. This argument assumes, of course, 
that our result is also relevant for the 0(N) model, but we are confident that this 
is actually the case. It would be nice if one could argue that the perfect classical 
action should then work equally well for nonabelian gauge theories. Unfortunately, 
already the large N limit of nonabelian gauge theories is so complicated that a 
similar analytic study is beyond our technical abilities. Finally, we like to comment 
once more on scaling versus asymptotic scaling. As our results show (and as various 
people have realized before) one should ask for scaling not for asymptotic scaling. 
In particular, the perfect action represents exact continuum physics although the 
bare lattice coupling does not exactly follow the perturbative /3-function. Still, for 
the perfect action asymptotic scaling sets in earlier than for the standard action. 

Although our real interest lies in QCD we believe that it may be helpful to 
gain further experience with fixed point actions for fermions in simpler settings like 
the Gross-Neveu model. For example, it would be interesting to investigate the 
fixed point action for a single set of staggered fermions (N = 2). Definitely, this 
is a numerical problem. Still, many results of the present paper are relevant for 
such an investigation. First, the locality of the fixed point action for free fermions 
has been optimized analytically resulting in c = 9 for the free parameter in the 
renormalization group transformation. This parameter should also be used when 
one looks for the perfect action at N = 2. Second, the small field fixed point action 
was derived at finite N. This is certainly the first step towards the full fixed point 
action. However, there is a lesson that we learnt: the small field fixed point action is 
not a perfect action. Hence, one should not rely on that approximation. Finally, the 
knowledge of the large N fixed point action may help to find a good parametrization 
for the perfect action at finite N. It is also important to explore fermionic fixed point 
actions in numerical simulations. In particular, it will be interesting to investigate 
how well standard numerical techniques like the Hybrid-Monte Carlo algorithm work 
in this case. 

It would be interesting to extend the large N analysis by systematically calcu- 
lating 1/N corrections to the fixed point action. This could be done not only for 
the Gross-Neveu model but also for 0(N) and CP(N) models. Of course, the real 
challenge is to develop perfect actions for the theory of strong interactions. 



29 



Acknowledgements 

We like to thank P. Hasenfratz, J. Jersak and F. Niedermayer for very interesting 
discussions, and J. Smit for a helpful remark. One of us (W.B.) likes to thank the 
Centro Brasileiro de Pesquisas Fisicas (Brazil) for warm hospitality. 

References 

[1] J. Kogut and K. Wilson, Phys. Rep. C12 (1974) 75. 

[2] K. Symanzik, Nucl. Phys. B226 (1983) 187, 205. 

[3] M. Liischer and P. Weisz, Commun. Math. Phys. 97 (1985) 59. 

[4] B. Sheikholeslami and R. Wohlert, Nucl. Phys. B259 (1985) 572. 

[5] W. Wetzel, Nucl. Phys. B255 (1985) 659. 

[6] P. Hasenfratz and F. Niedermayer, Nucl. Phys. B414 (1994) 785. 

[7] T. DeGrand, A. Hasenfratz, P. Hasenfratz, F. Niedermayer and U.-J. Wiese, 
work in progress. 

[8] D. Gross and A. Neveu, Phys. Rev. D10 (1974) 3235. 

[9] Y. Cohen, S. Elitzur and E. Rabinovici, Nucl. Phys. B220 [FS8] (1983) 102. 

[10] T. Jolicoeur, A. Morel and B. Petersson, Nucl. Phys. B274 (1986) 225. 

[11] L. Belanger, R. Lacaze, A. Morel, N. Attig, B. Petersson and M. Wolff, Nucl. 
Phys. B340 (1990) 245. 

[12] T. Kalkreuter, G. Mack, and M. Speh, Int. J. Mod. Phys. C3 (1992) 121. 

[13] W. Bietenholz and U.-J. Wiese, Nucl. Phys. B (Proc. Suppl.) 34 (1994) 516. 

[14] U.-J. Wiese, Phys. Lett. B315 (1993) 417. 

[15] J. Hirsch and S. Shenker, Phys. Rev. B27 (1983) 1736. 

[16] S. Coleman, Aspects of Symmetry, Cambridge University Press (1985). 



30 



